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We present results on semileptonic B decays obtained with the BABMl detector. The 
large data set accumulated at the PEP-II asymmetric-energy _B-Factory allows a new 
measurement technique, where the hadronic decay of one B meson is fully reconstructed 
and the semileptonic decay of the recoiling B meson is studied. Traditional analysis 
techniques of inclusive and exclusive B decays complement this approach with very 
high statistics data samples. These measurements play an important role in the tests 
of the description of CP violation in the Standard Model: The determinations of the 
Cabibbo-Kobayashi-Maskawa matrix elements \ Vab\ and \ Vub\ provide constraints on the 
unitarity of the CKM triangle. Furthermore, the experimental measurement of parame- 
ters of Heavy Quark Effective Theory test the consistency of the theoretical description 
of semileptonic B decays. 
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1. Introduction 

The principal motivation for the study of flavor physics is a comprehensive test of 
the Standard Model description of CP violation. Semileptonic B decays allow for 
a direct determination of \Vcb\ and \Vub\^ two elements of the Cabibbo-Kobayashi- 
Maskawa (CKM) quark mixing matrix. In the unitarity triangle, the precision of 
\Vcb\ affects constraints derived from kaon decays and the overall normalization, 
while the uncertainty in | Vub \ dominates the error of the length of the side opposite 
the angle (3. As this angle can be measured very cleanly in time-dependent CP 
asymmetries, the errors of \Vub\ must be model independent, well understood, and 
small before any discrepancies between sides and angles could be interpreted as new 

physics. This is not yet the case: the error in \Vub\ is larger than 10% and dominated 
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by theoretical uncertainties. 

In the theoretical description of semileptonic B decays, the large mass of the h 
quark plays a central role by implying special symmetries and a hard scale.^This is 
formulated by Heavy Quark Effective TheorjH (HQET) for exclusive decays and an 
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Operator Product Expansior!^ (OPE) for inclusive decays. Both provide systematic 
expansions of the (differential) decay rate in terms of Aqcd/ttii, and as{m{,). Here 
incalculable quantities are parametrized in terms of expectation values of hadronic 
matrix elements, which can be related to the shape (moments) of inclusive decay 
spectra. The large rate of Cabibbo-favored decays B Xctv allows for precise mea- 
surements of the relevant distributions and the determination of HQET parameters. 
This provides precision determinations of \Vcb\ and stringent quantitative tests of 
the consistency of the theory. 

The situation is different for Cabibbo-suppressed B Xutv decays: the large 
rate for B —f Xcdv decays constitutes a background that is about 50 times larger, 
overlapping in most of the phase space. Experimentally, selection criteria are applied 
to reduce this background, but can lead to problems in the theoretical description. 

The experimental approach to semileptonic decays can be separated into two 
classes: Exclusive decays reconstruct one signal decay mode, e.g., i?° D*^£~D. 
Even though the neutrino cannot be measured directly, this approach is relatively 
straightforward. The inclusive analysis of semileptonic decays is often based on the 
measurement of the charged lepton alone, or by combining hadronic final states X 
without disentangling specific resonances. 



2. The BABAR Detector 

The measurements presented here are based on data collected by the BABAR 
detectoJ'^ at the PEP-II asymmetric energy e+e~ collider near the T(45') resonance. 
Most of the analyses use an integrated luminosity of about 80fb~^, corresponding 
to about 89 million BB pairs. The T{AS) resonance is just above threshold for the 
decay into a pair of B mesons (either B^B~ or B'^B^), without any other fragmenta- 
tion particles. Furthermore, the two B mesons have a low and known momentum of 
p*g — 320MeV/c in the center-of-mass system (CMS)'', leading to spherically sym- 
metric decays. This is different for qq continuum processes (where q — u,d,s,c), 
which exhibit a more jet-like structure. This is exploited with event shape variables 
and neural networks. 

At a CMS energy of ^/s = 10.58 GeV, the T(4S') production cross-section 
amounts to about 1.1 nb. This corresponds to a rate of about 10 BB pairs/sec at 
an instantaneous luminosity of 10"^* cm^^ s~^. Given the total hadronic continuum 
cross section of ca. 3.5 nb, the resulting signal to background ratio is much higher 
than at hadronic colliders. The background from continuum processes is determined 
in dedicated "off-peak" runs, where the CMS energy is lowered to -^/s = 10.54 GeV. 

A five-layer silicon vertex tracker provides precision vertexing and low- 
momentum charged particle tracking, down to transverse momenta pj_ ^ 50 MeV. 
This is especially important for the reconstruction of D'^T^t decays, where 

the "slow" pion tt^ has very low energies. A 40-layer driftchamber surrounds the ver- 



'^AU variables measured in the CMS frame, e.g., p*g, are marked with a star. 
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tex detector and complements the momentum measurement. In addition, the dE/dx 
measurements are used in the identification of low-momentum electrons. The DIRC 
provides detection of internally reflected Cherenkov light used in charged hadron 
identification. The electromagnetic CsI(Tl) crystal calorimeter is the most impor- 
tant detector for electron identification (by means of the ratio E/p of the deposited 
energy E and the momentum p) . In addition, its measurements of neutral particles 
is crucial for the inclusive determination of the invariant mass mx in i? — > Xii? 
decays. The detector is surrounded by a superconducting coil (providing a magnetic 
field of 1.5 T) and its instrumented flux return, used in the identiflcation of muons. 

The boosted CMS at BABAR leads to a limited coverage of about 85% of the 
solid angle in the CMS. This is a notable disadvantage for the reconstruction of 
neutrinos from the missing momentum, as about 1 GeV of energy is missed per 
event (on average). 

3. Recoil Physics 

The very large luminosity at the BABAR detector allows for a new paradigm for 
the systematic study of semileptonic B decays. Traditionally'^ events are selected 
("tagged") by a high- momentum lepton, signaling the semileptonic decay of one of 
the B mesons and thereby reducing qq continuum events. 

At BABAR, an alternative event tagging technique has been developed: the 
hadronic decay of one B meson (Breco) is fully reconstructed and the semilep- 
tonic decay of the other B meson is identified by the presence of an electron or 
muon. This approach results in a low overall event selection efflciency, but allows 
for the determination of the momentum, charge, and flavor of the B mesons. It 
also provides a direct determination of the hadronic final state in B ^ XlD de- 
cays, as all particles in the recoil of the Breco candidate originate from the other B 
meson decaying semileptonically. This method also offers a promising way to study 
semileptonic B Xrvr decays. 

A very large sample of B mesons is reconstructed by selecting hadronic 
decays'^ Breco ~> DY^ ,D*Y^ , where the hadronic system consists of 

riiTT^ n2K^ n^K^ u^-k'^ , with ni + n2 < 5, nj, < 2, and rt4 < 2. The kinematic 
consistency of Breco candidates is checked with the beam energy-substituted mass 
jTiES = — p| and the energy difference /S.E = Eb — \/s/2, where ^/s is the 

total energy and (Eb,Pb) denotes the momentum four- vector of the Breco candi- 
date in the CMS. For each of the reconstructed B decay modes, the purity V is 
estimated as the signal fraction in events with tties > 5.27 GeV/ (see Fig. A 
priori, the purity of this sample is low, but improves substantially in conjunction 
with the requirement of a high momentum lepton in the recoil. By combining more 
than 300 modes, at least one Breco candidate is reconstructed in 0.3% (0.5%) of the 
B^B^ {B^B~) events. In events with more than one reconstructed Breco candidate, 

''Charge conjugation is implied throughout the text. 
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we select the decay mode with the highest purity. 




Fig. 1. The jties distributions for fully reconstructed hadronic B decays used in the event selection 
for the study of semileptonic B decays, a) With no requirement on the rest of the event, the purity 
amounts to 26%. b) With a p* > 1 GeV/c lepton, the purity improves to 67%. The arrows indicate 
the minimum ttses requirement used in the selection of signal events. 



4. B° D*+£-D 

While the decay mode D*^£~v has a large branching fraction, the measure- 

ments so far are not very consistent — recent results range from (4.59 ± 0.46)% to 
(6.09 ± 0.44)%. The very large luminosity at BABAR opens new possibilities for the 
precision determination of this decay. 

The theoretical description of i? ^ D^*^£D decays in terms of HQElEl predicts 
the differential decay rate schematically as 

dT{B^D^*)£D) r(u;2_ 1)1/2. 

-fi-\Vcb\ ■ { ,„2 1^3/2 -r■2/„,.^ (1) 



dw ~ ' \(w2-l)3/2.^2(^) 

where w = vb'Vd' — Ed* /mD* , J-'(^)(w) is the formfactor describing the hadroniza- 
tion into a meson, and /C is a known and constant factor. The Lorentz factor 
w of the c-quark in the 6-quark rest-frame takes values between w — \ ("zero- 
recoil" situation: the c-quark is at rest) and w = l.b (c-quark and £v leaving 
back-to-back). In the limit of toq — *■ oo the formfactors (w) are equal to the 
Isgur-Wise function.l^ Heavy quark symmetry provides the normalization constraint 
JF(^-)(1) = 1. At zero-recoil, the hght degrees of freedom (the spectator quark, the sea 
quarks and gluons) are not sensitive to the flavor change of the heavy quark. Because 
of the finite mass of the b and c quarks, small corrections need to be computed — 
this is done with phenomenological models or (currently quenched) lattice QCD 
calculations. 
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The experimental approach consists in the measurement of the differential rate 
dV/dw as a function of w, and the extrapolation of the data io w — 1 to obtain 
jF(l)|l/ct|. This measurement is preferentially done with B D*lv instead oi B ^ 
DiD decays: The decay rate is kinematically suppressed at w = 1 for both decays, 
but less so for B D*£v. By virtue of Luke's theorem^ there are no corrections 
at order 1/mf, or 1/toc for B D*£D, but they are present for B Dlv and 
thus increase the theoretical errors in this case. The background from high-mass Xc 
states constitutes a difficult experimental systematic problem. The decay B Dlv 
is even more affected by this, as the decay B D*£D here is a background process 
with a branching fraction that is abo ut t wo times larger. 

The event selection in this analysi^^ starts from a charged lepton (e or /z) with 
momentum p* > 1.2 GeV and a reconstructed D*+ — > D^tTs decay, where the D'^ is 
reconstructed in four modes: Z?" K^tt'^ , Kgir^n^ , K^tt^tt^tt^ , K^tt^tt'^ . The 
mass difference Sm = mjjo^^ — mjjo is used for the selection of !?*+ candidates 
and the determination of the combinatorial background. Given the very large data 
sample, it is possible to study and constrain most of the backgrounds directly in 
data: We study the uncorrelated background (where the lepton and D*+ originate 
from different B mesons) in control samples based on the opening angle between the 
D*~^ and the lepton (signal decays tend to be back-to-back). Continuum background 
is reduced with event shape variables, the remaining component is subtracted with 
off-resonance data. The determination of the most dangerous background from high- 
mass Xc states in data is described in the next section. The only component taken 
from Monte Carlo (MC) simulations is the correlated background, where B 
D*+X,X -> Ye. 

At the T(45), the known momentum magnitude of the B mesons provides sen- 
sitivity to the missing mass in the signal decay 5° D*~^e~D from the observed 
particles D*~^ and £. Assuming that the only unmeasured particle is a massless 
neutrino, the angle between the momenta of the B'^ and the combined D*£ is 

COSOBM'i- j ■ l^J 

2\pb\\pd'i\ 

This quantity will lie in the physical range for a signal decay (modulo resolution 
effects; signal events are required to have | cos ^b.D'^I < 1-2). Decays with additional 
missing particles will lead to a tail at low values, illustrated in Fig. [21 for B'^ — ^ 
D*'^e~i/ decays, where the background from D**iD exhibits a long tail. This is 
also visible for signal decays, due to missed photons from bremsstrahlung. For the 
extrapolation to zero-recoil, the decay rate must be measured differentially in w = 
(m^o — q^)/{2mBomD'), where q"^ — {pgo — po')^ ■ The direction of the B^ 

momentum is obtained from eq.Elup to an azimuthal ambiguity about the direction 
of the D*e pair: an unbiased estimator of w with a resolution of a{w) ~ 0.02 is 
calculated from the average of the two solutions corresponding to minimal and 
maximal angle between the and D*+ mesons. 

The signal yield as a function of w is determined with a fit based on a quadratic 



6 Urs Langenegger 



formfactor parametrizatioii^ (see Fig. Based on ca. 57000 signal events, we 
obtain |Vch| = (37.27 ± 0.26stat ± lASsysttHltheo^ ^ dominant system- 

atic errors of this result are due to tracking and vertexing (but not the slow pion 
efficiency), branching fractions and B lifetimes, and /oo = B{T{AS) — > B'^B^). 
The theoretical error in the lattice calculation for i^*(l) = O.OlSlgQgg is balanced 
among statistical, fitting, matching, spacing, mass, and quenching components-l^ 




COsQgy^^l 1 1.05 1.1 1.15 1.2 1.25 1.3 1.35 1.4 1.45 1.5 

w 



Fig. 2. a) Distribution of cosff^n ]jtf for the decay — » D*+e~i7e,D'' K~tt^, where the 
points show the data and the filled histograms represent the MC components. The MC normal- 
ization is determined in a fit. The bottom plot shows the fractional deviation of the data from the 
fit result, b) Comparison of the w distribution for f?" — > D'+e"!/^ in data and the result of the 
fit. The fit residuals are shown in the bottom plot. 



The branching fraction B{B D*lv) = (4.68 ± O.OSstat ± 0.29^yst)% is deter- 
mined by integrating the differential w distribution (cf. cq. ^ thus reducing the 
uncertainties from formfactor parametrizations. These results are somewhat lower 
than other measurements (especially the recent one by CLECJ^D. As all measure- 
ments are systematically limited, a complementary approach, e.g., in the recoil of 
a fully reconstructed B candidate, will be a very interesting result addressing the 
largest errors. The theoretical uncertainties associated with the extrapolation to 
w = 1 might be much reduced with a sample of the order of 10^ signal decays. 

5. Inclusive Cabibbo-favored Decays 

In the description of exclusive B^ — > D*^t~v decays, the HQET relies on both 
rrib and rric being large. For inclusive decays, this requirement can be relaxed to 
rrih ^ ^QCD, expressing the separation between the very short time scale rele- 
vant for the weak 6-quark decay and the long time scale for the hadronization of 
the hadronic remnant. The OPE^can be combined with HQET to calculate, e.g., 
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the total semileptonic width T si in a power series in Aqcd I'^h and a perturbative 
expansion in OLs{vn\y\ Incalculable quantities are parametrized in terms of nonper- 
turbative hadronic matrix elements. At lowest order, the OPE expression reduces 
to the parton model. There are no power corrections at order v1qcd/"16 in the total 
rate. The leading corrections are parametrized with 

• yl = TOB — Wf, + (Ai + 3A2)/2mb, the energy of the light degrees of freedom. 
To first order, A is the mass difference of the h quark and the B meson. 

• Ai is the negative kinetic energy squared of the h quark in the B meson. 

• A2 describes the chromomagnetic coupling of the h quark spin to the light 
degrees of freedom. 

Note that these parameters are scheme and order dependent. At higher orders 
many more parameters (pi, p2,^, • ■ •'?4, etc.) enter, none of which are currently 
well known. They constitute a large fraction of the theoretical errors. 

These parameters are not restricted to the description of the total semileptonic 
rate, they also appear in the calculation of differential semileptonic B decay spectra 
and in other B decays: A can be related to the mean photon energy of the decay 
h 57, and A2 can be determined from the mass difference of B* and B mesons. 

The parameters can be related to the shape of decay spectra in semileptonic 
B XiD decays, such as the lepton energy spectrum Ei or the invariant hadronic 
mass spectrum mx- The inclusive description of semileptonic B ^ XiD decays 
does not distinguish between specific hadronic final states X, e.g., the D and D* 
mesons. To compare, e.g., the theoretical expectation for the mx distribution with 
the measured spectrum, it is therefore necessary to resort to observables smearing 
the differential spectrum. A simple and sensitive possibility is given by moments of 
various order. As different moments have different dependencies on the parameters, 
a simultaneous fit to several moments provides for an experimental determination of 
the nonpcrturbative parameters. The large branching fraction for B Xii> allows 
for a very precise determination of these parameters. We can therefore shift a large 
fraction of the theoretical errors into (smaller) experimental errors, which are more 
amenable to a proper statistical interpretation. This is essential for a quantitative 
understanding of the errors of the extracted CKM parameters. 

5.1. Inclusive Semileptonic Branching Fraction 

The model-independent measurement of the total inclusive semileptonic branching 
fraction B{B Xe~v) was pioneered by the ARGUS collaboration!^ It has a small 
model-dependence in the sense that no assumptions on the shape of the primary 
electron spectrum from B — > Xe~iy decays are necessary. In this analyis,^^we select 
events by the presence of a high-momentum tag electron (1.4 < p* < 2.3 GeV). In 
the rest of the event, signal electrons are identified and grouped into two separate 
classes, depending on whether they have opposite charge ("unlike sign sample" : the 
two electrons are either from B X^e^v, B — s- X^e^v or from B — s- XcE^v, Xc 
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Yci^v) or the same charge ("like sign sample": B Xce^v, B ^ B ^ XcC^v and 
B XcC^v, B Xc, Xc Ye^v) as the tag electron. 

In the unlike-sign class, the background B Xce~v,Xc Ye^v can be 
strongly reduced by exploiting the fact that electron pairs from the same B are 
preferentially back-to-back, while there is no correlation in the opening angle distri- 
bution a{etag,e signal) for the case of the pair coming from two different B mesons. 
This is illustrated in Fig. where the opening angle distribution for signal elec- 
trons with 0.7 < p* < 0.8 GcV shows a flat signal component and a background 
contribution rising toward low values. The shape information is taken from MC 
simulation, and the data distribution is fitted to determine the background con- 
tribution (shaded) in the signal region. The requirement of a small opening angle 
also removes most of the heavy pair background (from J/ijj — > e^e^\ the remain- 
ing pair background from photon conversions and Dalitz tt*^ e^e^7 decays is 
reconstructed and removed explicitly). The effect of B^B^ mixing can be unfolded, 
as 



1 dN-^z!^ C?iVp^2mary 



(l-x) + ^^-X (3) 



where X = fooXd = 0.087 (here Xd = 0.174±0.009 is the B°B° mixing parameteJ^ 
and foo = 0.50 has been assumed) and the efficiency £a(p*) of the opening angle 
requirement. 




cosa P [GeV/c] 

Fig. 3. a) Cosine of the opening angle between the signal electron with 0.7 < p* < 0.8GeV/c and 
the tag electron in the unlike-sign sample. The shaded area represents the background electrons, the 
vertical line illustrates the minimum requirement on the opening angle, b) Momentum distribution 
of electrons from B — > Xe~i' decays, after efficiency and bremsstrahlung corrections. 



The integration of the spectrum over the range 0.6 < p* < 2.5 GeV (see Fig. Ob) 
yields N{B Xe~v) — 25070±410stat signal events for an integrated luminosity of 
about 4 fb^^. Only small corrections need to be applied: Bremsstrahlung corrections 
(2.2%), geometric acceptance (16%), event selection bias (2%) and the extrapolation 
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(6.1%) to p* = 0. From these numbers and the overall normalization from the 
number of tag electrons we determine B{B Xe'v) = (10.87±0.18±0.30)%. The 
dominant systematic error arise from electron identification plus tracking and from 
semileptonic decays of upper-vertex charm particles, especially affected by the poor 
knowledge of B{Ds — > 07r). Extending the measurement range to lower values of p* 
does not improve the error, as backgrounds (pair background, cascade decays) with 
large uncertainties grow very large. 

From the total branching fraction, \Vcb\ can be extracted, e.g., in the IS 
expansion.^! The errors in this approach are dominated by theoretical uncertain- 
ties. In the next section an alternative method is described that takes into account 
more information and their correlations. There, the differential measurement of the 
lepton energy spectrum and its moments will contribute substantially. 

The measurement of i? — > Xe^v has also been done in the recoil of a Breco 
candidate,^! albeit with much less statistics. This allows a comparison of electron 
energy spectra from B'^ and B^ decays and provides one way to study effects of 
quark-hadron duality violation and other nonperturbative effects like weak annihi- 
lation and Pauli interference. 



5.2. Hadronic Mass Moments 

In B XiD decays, the invariant mass mx distribution is the most sensitive probe 
to physics beyond the parton model and hence to the nonperturbative parameters 
A and Ai. The lepton energy distribution also has sensitivity, but at a reduced level. 
The experimental feasibility matches this situation: the lepton energy distribution 
can be obtained with high precision and resolution from a measurement of the elec- 
tron only, a relatively straightforward task. On the other hand, the reconstruction 
of the complete hadronic final state i s a much more involved procedure. 

The event selection in this analysiff^ requires a Breco candidate and an identified 
lepton with p* > 900 MeV/c and charge consistent for a primary B decay. The charge 
imbalance of the event is required to be not larger than one. These criteria lead to 
a data sample of about 7100 events. 

All remaining charged tracks and neutral showers that are not part of the Breco 
candidate are combined into the hadronic system X. A neutrino candidate is recon- 
structed from the missing four-momentum pmiss = Pt{4S) ~ Px — PB^^^oi where all 
momenta are measured in the laboratory frame. Consistency of the measured Pmiss 
with the neutrino hypothesis is enforced with the requirements Emiss > 0.5 GeV, 
IPmissI > 0.5 GeV, and \Emiss — IPmissW < 0.5 GeV. The determination of the mass 
of the hadronic system is improved by a kinematic fit that imposes four-momentum 
conservation, the equality of the masses of the two B mesons, and forces pi = 0. 
The resulting mx resolution is 350 MeV/c^ on average. MC simulated event samples 
are used to calibrate the absolute mass scale, determine efficiencies, and estimate 
backgrounds. This mass scale calibration allows the direct determination of the 
moments of the mx distribution without recourse to simulated decay spectra. 
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The resulting moments of the hadronic mass-squared distribution are shown as 
a function of the threshold lepton momentum p^j„ in Fig. A substantial rise 
of the moments toward lower momentum is visible, due to the enhanced contribu- 
tions of high-mass charm states (phase-space suppressed at higher J3^i„). The main 
contributions to the systematic error are from the detector response simulation and 
from semileptonic decays of upper-vertex charm particles. The uncertainty from the 
modeling of the Xc state is negligible compared to the other systematic errors. 
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Fig. 4. (a) Measured hadronic mass moments for different lepton threshold momenta p^j„. The 
errors of the individual BABAR measurements are highly correlated. For comparison, the measure- 
ments by the DELPHI and CLEO collaborations are also shown. The solid curve is a fit to the 
BABAR data; the dashed curve is the OPE prediction based on the CLEO result combined with 
information from the decay b — > 57. (b) Constraints on the b quark mass and \Vcb\ from this 
measurement, and the fit to the combined hadron moments and lepton moments, respectively. 



Accounting for all correlations between the moments at different p^j„, we ob- 
tain A = 0.53 ± 0.09 GeV and Ai = -0.36 ± 0.09 GeV^ in the MS regularization 
schemc^^The errors given do not include uncertainties due to terms at 0{l/m^). 
For comparisOTi, Fig. 0^ also shows the result of the hadronic massmeasurement 
of DELPHI,'^ fully consistent with this result. The CLEO resullP2 of the first 
hadronic mass moment at = 1.5 GeV is also consistent, but in combination 

with the mean photon energy from b s-^l^ shows a different dependence 
(taking into a ccou nt the bias from the limited photon energy range, the agreement 
is much betteJSHl). 

A fit to all hadronic moments from BABAR is performed in the IS* scheme,!^ as 
this scheme exhibits better convergence of the perturbative series than other alter- 
natives. The results are = 4.638 ± 0.094e:Ep ± 0.062dimesLM ± 0.0651/^1 GeV 
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and Ai = — O.26±O.O6ea;p±O.O4dim0BLM±O.O4i/„j| GeV^. The fit also utilizes the 
semileptonic width Tsi = (4.37 ±0.18) x 10""MeV (determined from BABylRdata) 
to determine \Vcb\ = (42.10 ± IMexp ± 0.52dim(BBLM ± 0.50i/„,3^) x 10-^. 

The consistency of the OPE is tested by combining the measurement of BABAR 
with the four lepton energy moments measured by the CLEO collaboration-^^ and 
the DELPHI collaboration. ^0 In Fig. ^p, the fit results are shown separately for 
hadron mass and lepton energy moments. The Ax^ = 1 contours of hadronic mass 
and lepton energy moments do not overlap. Note that the theoretical errors here 
do not include uncertainties due to terms at 0(l/m^). In the future, the exten- 
sion of these measurement to include more high-precision observables will allow a 
quantification as to what precision can be expected realistically from the OPE. 

6. Inclusive Cabibbo-suppressed Decays 

In the measurement oi B X^tv decays, the large background from B Xc^v de- 
cays has to be reduced by restricting the phase space in the analyses. One possibility 
is to measure the lepton energy spectrum at the "endpoint" , beyond the kinematic 
cutoff for B Xctv decays, at p* > 2.3 GeV. A disadvantage of this approach 
is that only about 10% of all charmless semileptonic decays are measured. This 
leads to a significant extrapolation with corresponding uncertainties. The model- 
dependence of this error can be reduced with information on the movement of the 
b quar k in side the B meson obtained from the photon energy spectrum in 6 — > 37 
decays .12^1 

Alternative methods have been proposed, such as the invariant mass of the 
hadronic system X in B ^ XiP decays.!^ Here 50-80% of all B X^tv decays 
are measured, depending which cut mx < m'^* is used for the determination of the 
signal yield. As in the case of the endpoint spectrum, there is a dependence on the 
light-cone distribution function ("shape function") of the b quark in the B meson, 
describing the Fermi motion of the heavy quark in the meson. Another process sen- 
sitive to this shape function is the rare decay b 57, with a branching fractiorP2 
B{b sj) = 3.21 ± 0.53 x lO""*. This rate is even smaller than B{B XuiP), thus 
limiting the precision of the determination of the shape function parametrization 
and parameters. If the HQET parameters determined in Cabibbo-favored semilep- 
tonic B decays can be used consistently in the shape function description, the 
determination of \Vub\ will benefit from the high precision measurements described 
in the previous section. 

The total rate can be translated to \Vub\ with an error of about 5% from uncer- 
tainties of higher orders in the perturbative expansion and the uncertainty of the b 
quark mass .'12'^ 

6.1. Endpoint Spectrum 

In this analysis ,12^ the event selection is based on a high-momentum electron 
{p* > 2.0 GeV) and the signature of a neutrino. Specifically, we require for the 
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missing momentum Pmiss > 1 GeV to be pointing into the main detector accep- 
tance —0.9 < cos6'*„jgg < 0.8 and in the hemisphere opposite to the electron. At 
low momenta, the electron sample is dominated by electrons from B XcHu de- 
cays over a continuum background. This latter component dominates the spectrum 
at high momenta. We fit the continuum background with a 4*'' degree Cheby- 
shev polynomial in the off-resonance data and the high-momentum range in the 
on-resonance data. After subtraction, the signal events are visible in the range 
2.3 < p* < 2.6 GeV and illustrated with the solid (red) histogram in Fig. ^p. The 
restriction to this momentum range yields a total of 1696 ± 133 signal events with 
a signal to background ratio S/B ~ 0.25. Extending the momentum range to 
lower values decreases S/B due to more background from B X^iv decays, lead- 
ing to substantially higher uncertainties from the modeling of -B — > XcH-v decays. 
On the other hand, the extrapolation is decreased, resulting in a smaller theoret- 
ical error. In the endpoint range, the partial branching fraction is determined to 
be AK(i? XuiD) = (0.152 ± O.OU^t^t ± O.OOU^y^*) x lO'^, where the dominant 
errors arise from the uncertainties in the continuum subtraction, the motion of the 
B meson in the T{AS) rest-frame, and the selection efficiency. 

From this result, the extrapolation to the total semileptonic charmless branching 
fraction is done as in the CLEO analysis. Here the shape function parameters are 
determined by a fit to the b s-f photon energy spectrum. The result is B{B 
Xutv) = (2.05 ± 0.27ea;p ± 0-46/„) x 10^'^, where the errors are now grouped into 
a first part containing the statistical and systematic uncertainty from the endpoint 
measurement and a second part describing the uncertainties from the extrapolation. 
This yields iKhl = (4.43 ± 0.29e^p ± 0.50/„ ± 0.35^^ ± 0.25r) x 10"^. 



6.2. Hadronic Mass Spectrum 

In this analysiJ2^ we use the invariant mass mx of the hadronic system to sepa- 
rate B Xutv decays from the dominant B Xctv background in events tagged 
by the fully reconstructed hadronic decay of a Breco candidate. This method of- 
fers a substantially larger signal acceptance than the endpoint measurement. The 
hadronic system X in the decay B X£D is reconstructed from charged tracks 
and energy depositions in the calorimeter not associated with the Breco candidate 
or the identified lepton. We require exactly one charged lepton with p* > 1 GeV/c, 
charge conservation {Qx + Qe + Qb^^co = 0)i "^miss < 0.5 GeV^. We reduce the 
B'^ — > D*^t^T' background with a partial reconstruction of the decay (using the 7r+ 
from the _D*+ D'^'^t decay and the lepton). Furthermore, we veto events with 
charged or neutral kaons in the recoil B. 

In order to reduce experimental systematic errors (in particular lepton identi- 
fication), we determine the ratio of branching fractions Ru from N^, the observed 
number of B Xu£i> candidates with mx < 1.55 GeV/c^, and N^i = 29982 ± 233, 
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Electron Momentum (GeV/c) Electron Momentum (GeV/c) 

Fig. 5. a) Electron momentum distribution in the T(4S) rest-frame for on-resonance and off- 
resonance data, respectively, b) Electron momentum distribution after continuum subtraction, 
with signal and background MC distributions, c) The differential branching fraction as a function 
of the electron momentum, after efficiency and bremsstrahlung corrections. The data are compared 
to the prediction of the ISGW2 model (assuming a total inclusive branching fraction of 10^'' for 
B — » Xuiv decays with mass mx < 1.5GeV/c^). 



the number of events with at least one charged lepton: 

^ _ B{B ^ X^iD) _ N^/{e^,ie-^,) efef,,. 

Here e"^; — 0.342 ± 0.006sta( is the efficiency for selecting B X^iD decays once a 
B Xiv candidate has been identified, ejjj^ = 0.733 ± 0.009stat is the fraction of 
signal events with the reconstructed mx < 1-55 GeV/c^; sf/ef — 0.887 ± O.OOSstat 
corrects for the difference in the efficiency of the lepton momentum cut for B — > XiD 
and B XuiD decays, and sfeco/^reco = 1-00 ^ 0-03stat accounts for a possible 
efficiency difference in the Breco reconstruction in events with B —f XiD and B — > 
Xuii' decays. 

We extract Nu from the mx distribution by a fit to the sum of three contri- 
butions: signal, background Nc from B — > Xcii', and a background of < 1% from 
other sources. In each bin of the mx distribution, the combinatorial Breco back- 
ground for Toes > 5.27 is subtracted on the basis of a fit to the toes distribution. 
Fig. EJi shows the fitted mx distribution. To minimize the model dependence, the 
first bin is extended to mx < 1.55GcV/c^. We find 175 ± 21 signal events and 
90 ± 5 background events in the region mx < 1.55 GeV. From this wc determine 
Ru = (2.06 ± 0.25stat ± 0.23syst ± OMtheo) x 10"^ The dominant detector sys- 
tematic errors are due to the uncertainty in photon detection and combinatorial 
background subtraction. The efficiencies e"g; and e"!^^ are sensitive to the model- 
ing of the B Xu(i> decays.l^ We assess the theoretical uncertainties by vary- 
ing the nonperturbative parameters within their errors, A — 0.48 ± 0.12 GeV and 
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Ai = —0.30 ± 0.11 GeV^, obtained from the results in Ref. EJby removing terms 
proportional to 1/m^ and aj. from the relation between the measured observables 
and A and Ai. Here we assume that the parameters of the shape function are given 
by the HQET parameters. 
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Fig. 6. The mx distribution for B — > Xtv candidates: a) data (points) and fit components, and 
b) data and signal MC after subtraction of the b — > ctv and the "other" backgrounds. 



Combining the ratio Ru with the measured inclusive semileptonic branching 
fractionfl^i we obtain SfB XJD) = (2.24 ± 0.27 ± 0.26 ± 0.39) x lO'^. With 
the average B lifetimtEII we obtain \Vub\ = (4.62 ± 0.28stat ± 0.27syst ± 0.40t/ieo ± 
0.26r) X 10~'^, where the last error is the uncertainty in the extraction of \Vub \ from 
the total decay rate. No error is assigned to the assumption of parton-hadron duality. 
This result is consistent with previous inclusive measurements, but has a smaller 
systematic error, primarily due to larger acceptance and higher sample purity. 

7. Exclusive Cabibbo-suppressed Decays 

The reconstruction of exclusive charmless semileptonic B Xuli' decays, where 
Xu = TT, p, w, . . ., is challenging due to large backgrounds and the missing neutrino. 
The theoretical description is on less solid foundations than in the inclusive case. 
The determination of \Vub\ in this case requires the computation of formfactors, 
parametrizing the behavior of the hadronic current in the B meson decay. Due to 
the light mass of the hadronic final state, HQET is of much less help here compared 
to -B — > D*£v, and most of the calculations are based on phenomenological models. 
Currently, lattice QCD calculations are still based on the quenched approximation, 
with uncertainties of the order of 15-20%. The decay B — > iriv is more amenable 
to these calculations than the decay B — > plv, where the hadronic final state is a 
broad resonance. The computation of formfactors is only possible in the region of 
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low momentum Ptt < 1 GeV. This corresponds to high > 17 GeV^, where the rate 
is kinematicaUy suppressed. In addition, this is precisely the kinematic region where 
the current experimental methods are limited by very high backgrounds. With the 
high luminosities achievable at the B factories, the measurement of exclusive decays 
in the recoil of a Breco candidate offers large advantages. 

7.1. B° p-e+9 

This analysiJ^ aims for a h igh neutrino reconstruction efficiency and is similar to a 
previous CLEO analysis .^3 Events are selected with a high- momentum electron and 
divided into two samples based on the electron momentum: The high-momentum 
high-Ee with 2.3 < < 2.7 GeV and low-Ee 2.0 < E^ < 2.3 GeV. The high-Ee 
sample is primarily used for the determination of the signal, while the low-Ee sample 
serves for the measurement of the B —^ Xc(.D background. The analysis is optimized 
for the measurement of B — + pe+i^, but also selects B — > Tre^v and B^ — > uje^v to 
better control cross- feed background contributions. 

Continuum background is suppressed with a neural network of 14 variables. As 
the hadronic final state in this exclusive analysis is much less constrained than in 
B — > D*£D, the neutrino reconstruction carries more weight. The direction of the 
missing momentum is required to point into the main detector {\ cos 9miss\ < 0.9) 
to reject events with substantial energy loss along the beam axis. The angle a 
between the reconstructed missing momentum and the inferred neutrino momentum 
is required to be small (cos a > 0.8) and (as in the decay B — > D*£D) the angle 
between the B meson momentum and the combined eh momentum (where h = 
Tr,p,u!) is required to lie in the physical region {\cos0B,eh\ < l-l)- 

The signal yield is extracted by a binned maximum likelihood fit to high-Ee 
and low-Ef, samples in two variables: AE* — E^ad + E^ + pp — Ebeam (where pp is 
obtained from the missing momentum) and the mass of the hadronic system nihad = 
m^7r(7r)- The shape of the continuum distributions is taken from off-peak data, the 
remaining shapes are from MC simulations. The fit incorporates isospin and quark 
model relations. Fig. [Ji illustrates the m7r,r variable in the high-Ee sample. For an 
integrated luminosity of £ = 50fb~^, a signal yield of 5' = 505 ± 63stat events is 
obtained, leading to a branching fraction of S(_B" p^e^v) = (3.29 ± 0A2stat ± 
0A7syst ± 0.60t;ieo) X 10^*. The CKM matrix element \Vub\ is determined from the 
relation - \f^^^^§^ = (3.64 ± 0.22,^,* ± 0.25,^,*!° j^^^^J x IQ-^. The 
dominating systematic error are from the modeling of resonant and nonresonant 
B Xutv decays, the tracking efficiency, and the fit method. The theoretical 
error is determined as half of the full spread of all theoretical uncertainties in the 
formfactor calculations A measurement of the dependence will help to reject 
models not describing B — > pei> decays, but will not help in reducing the inherent 
model dependence of the error. Here, only unquenched lattice QCD calculation of 
the formfactors and the experimental measurement in the same kinematic region 
will advance the field. 



16 Urs Langenegger 



data events exclusive 




0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 ,2 
Mhad [GeV/c ] 




-0.3 -0.2 -0.1 -0 0.1 0.2 0.3 0.4 

m,(GeV/c=) 



Fig. 7. Illustration of tlie signal to background ratio for exclusive charmless semileptonic decays, 
a) Continuum-subtracted fit projections for m^o^- for the p~e'^v channel in the high-Ee 

electron-energy region. The contributions are the direct and crossfeed components of the signal 
(unhatched region, above and below the dashed line, respectively), the background from b — > ueu 
(double-hatched region), and b — > ceu and other backgrounds (single-hatched region), b) The mx 
projection for B~ — > ir^i~v measured on the recoil of a fully reconstructed hadronic B decay. 



7.2. Exclusive Decays on the Recoil 

This analysii^ is a combination of the high-purity event tagging based on a fuUy 
reconstructed hadronic B decay and the exclusive reconstruction of signal decay in 
the recoil. This approach results in a very low overall signal efficiency of the order 
of 0.1%, but allows the measurement of i? — > -Kiv over the entire kinematic range. 
This type of measurement will become a prime method for the determination of \ Vub\ 
from exclusive semileptonic decays, as the traditional approach mentioned in the 
previous section is affected by very low S/ B problems, especially in the range where 
lattice QCD can provide model-independent formfactor calculations. Because the 
statistical yield of the method is low, an integrated luminosity of about C ~ 500 fb^^ 
is needed before the method will provide a better measurement of _B ^ T:li> than 
the traditional approach. 

While all exclusive decays can be measured in this analysis paradigm, only 
the measurement B iriv shall be described here (for B^ Tr'^£~i>) as it 
is the most promising channel for lattice QCD calculations. After the require- 
ment of a fully reconstructed hadronic B decay for tagging purposes, events with 
semileptonic B decays are selected by the requirement of a muon or electron with 
pe > 1.0 GeV. Cabibbo-favored B Xctv decays are rejected by requirements 
on the missing mass rn^^^g < 0.4 GeV^, the invariant mass of the tt" candidate 
and the requirement that no other hadronic charged track be detected. The num- 
ber of signal S = 7.0 ± 2.6 events is determined from a fit to the tties distri- 
bution of selected events, corrected for background B = 0.2 ± 0.2 (determined 
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from MC simulations) and selection efficiency e — 0.42 ± 0.04 for all requirements 
after the Breco and lepton candidate selections. As in the case of the inclusive 
B — > Xuiiy analysis, the signal yield is normalized to the number of events with a 
charged lepton in the recoil of the Breco candidate. The result of this analysis yields 
B{B- Tr°e-i})/B{B XiD) = (0.76 ± OM.tat ± O.ll^j;^*) x lO'^. This result 
is statistics limited; the largest systematic errors are uncertainties of the tties fits, 
a possible selection bias for charmless semileptonic B decays compared to general 
semileptonic B decays, and the measurement of photons in the calorimeter. 

8. Conclusions 

In the last few years, the study of semileptonic B decays has offered many new 
perspectives. Theoretical uncertainties are parametrized in terms of experimental 
observables, substantially reducing the model-dependent component in the total 
error on \Vcb\- The large luminosity at experiments like BABAR at the T{4:S) res- 
onance allows to select events by means of fully reconstructed hadronic B decays 
and studying the semileptonic decay of the other B meson. This opens the precise 
study of spectral moments in semileptonic B decays and therefore the precision 
determination of \Vcb\ and the underlying fundamental parameters of the theory. 
It also allows for better constraints in the study of _B ^ Xuii^ decays and leads 
to improved determinations of \Vub\- In the future, the measurement of exclusive 
charmless semileptonic B decays in combination with unquenched lattice QCD cal- 
culations will challenge the precision of inclusive \Vub\ determinations. In the far 
future, leptonic B decays may provide yet another way to \Vub\- 
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